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Abstract. Let {px)\£A be a holomorphic family of representations of a surface group 7ri(S') 
into PSL(2,C), where 5 is a topological (possibly punctured) surface with negative Euler 
characteristic. Given a structure of Riemann surface of finite type on S we construct a bi- 
furcation current on the parameter space A, that is a (1,1) positive closed current attached 
to the bifurcations of the family. It is defined as the dd'^ of the Lyapunov exponent of the 
representation with respect to the Brownian motion on the Riemann surface 5", endowed 
with its Poincare metric. We show that this bifurcation current describes the asymptotic 
distribution of various codimension 1 phenomena in A. For instance, the random hypersur- 
faces of A defined by the condition that a random closed geodesic on S is mapped under 
Pa to a parabolic element or the identity are asymptotically equidistributed with respect to 
the bifurcation current. The proofs are based on our previous work [DDT] , and on a careful 
control of a discretization procedure of the Brownian motion. 



Introduction 

Let G be a finitely generated group and (pa)agA be a holomorphic family of representations 
of G into the Mobius group PSL(2,C). Viewing these representations as dynamical systems 
on P^, we can divide the parameter space A into a stability locus, where the action of G is 
locally quasi-conformally structurally stable, and its complement, the bifurcation locus. 

In a previous article [DDlj, under natural assumptions on the family (/Oa)agA) we showed 
how to design bifurcation currents on A. These are positive closed currents aiming at an 
understanding of the measurable and complex analytic structure of the bifurcation locus. This 
construction relies on probabilistic ideas, let us briefly recall some of its key steps. We start by 
choosing a probability measure ^ on G, which gives rise to a random walk on G. An example 
is the nearest neighbor random walk on the Cayley graph associated to some symmetric set 
of generators. Fix any matrix norm ||-|| on PSL(2, C). Then the theory of random products of 
matrices implies that for any non-elementary representation p, for almost every sample path 
[gn - ■ ■ gi)n>i of the random walk, ^\og\\p{gn ■ ■ ■ gi)\\ converges to the Lyapunov exponent 
Xnip)- It is easily seen that the Lyapunov exponent function A Xfi{p\) is plurisubharmonic 
(psh) on A, and we define the bifurcation current by the formula Tbif,^ = dd'^Xfj.iPx)- 

It turns out that the support of Tj^a^fj, is precisely equal to the bifurcation locus. Also, Tbif,^ 
describes the asymptotic distribution of certain dynamically defined sequences of subvarieties 
in A. Indeed, for t G C, and g G G, let Z{g, t) be the hypersurface in A defined by 

Z(5,t) = {AGA, ti\px{g))=t} 

(tr stands for the trace). The emblematic value is t = 4, corresponding to those parameters 
where px{g) becomes accidentally parabolic or equal to the identity, which are known to be 
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dense in the bifurcation locus. Then the following equidistribution statement holds: for almost 
every sample path of the random walk, the sequence of divisors ^Z{gn ■ ■ ■ gi,t) converges to 
Tbif,/x when n — )• oo. 

We note that analogous results were previously known in the context of holomorphic fam- 
ilies of rational mappings on P^. An essential motivation in |DDlj was actually to fill the 
corresponding line in the famous Sullivan dictionary between rational iteration and Kleinian 
group theory. The reader is referred to |Duj| for an overview on these issues. 

In this paper we specialize to the case where G = vri(S') is the fundamental group of 
an orientable surface of genus g with n punctures, such that xiS) = 2 — 2g — n < 0. As 
before we consider a holomorphic family {px)x(=A of non-elementary representations of G into 
PSL(2, C) (we can equally deal with representations modulo conjugacy) sending simple loops 
about the punctures to parabolic transformations. Such representations will be referred to 
as parabolic. For the sake of this introduction, let us assume without further notice that 
these representations are faithful at generic parameters and not conjugate to each other in 
PSL(2,C). Our main purpose is to construct a bifurcation current on A associated to the 
choice of a conformal structure on S. For this, let X be a structure of Riemann surface of 
finite type on S (that is, X is biholomorphic to a compact Riemann surface minus n points) 
and endow X with its unique hyperbolic metric of curvature —1. Recall that it is of finite 
volume. This data in turn gives birth to a wealth of mathematical objects: 

- a dynamical system: the geodesic flow on the unit tangent bundle T^X; 

- a stochastic process: Brownian motion on X; 

- an isomorphism between G and a lattice F in PSL(2, R) induced by the uniformization 
of X 

We will demonstrate that our bifurcation current fits naturally with all these structures. Let 
us already point out that as in the previous case, this bifurcation current will be defined by 
taking the dd'^ of a certain Lyapunov exponent function. 

A family of representations of G (modulo conjugacy) naturally associated to the Riemann 
surface structure is obtained by considering the holonomies of holomorphic projective struc- 
tures on X. The parameter space for this family of representations is biholomorphic to 
(j^3g-3+n ^^^^^ jg ^j^g ambient space for the so-called Bers embedding of Teichmiiller space. As 
such it plays an important role in Teichmiiller theory (see Dumas |Dumj for a nice introduc- 
tion to this topic). Notice that the much studied cases (5 = 0, n = 4) and (5 = l,n = 1), 
corresponding to one-dimensional Teichmiiller spaces, are associated to punctured surface 
groups. This is a motivation for allowing for non-compact surfaces, which, as we will see, are 
a source of technical difficulty. 

We see that the space of projective structures carries a natural bifurcation current (resp. 
Lyapunov exponent function) . The Sullivan dictionary predicts that the structure of this space 
should resemble that of the space of polynomials in holomorphic dynamics. In a companion 
paper jDD2j we introduce an invariant, the degree of a projective structure, and show that it 
is closely related to the Lyapunov exponent of its holonomy. This enables us to unveil a new 
analogy between spaces of projective structures and spaces of polynomials. 

Let us discuss more precisely the contents of the present paper. As before, the definition 
of the bifurcation current depends on the construction of a Lyapunov exponent function 
on A, that we now describe. Mark a point * in X, which will serve as a base point for the 
fundamental group, and equipp X with its Poincare metric. For a Brownian path uj : [0, 00) — )• 
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X starting at a;(0) = we do the following construction: for t > 0, we append w|[o,t] with 
a shortest path returning to to get a closed loop ujt, hence an element [oot] G 7ri(X, -we 
disregard any possible ambiguity in the choice of the returning path. 
We then have the following result 

Proposition A. Let X he a hyperbolic Riemann surface of finite type and p : G ^ PSL(2,C) 
be a non- elementary parabolic representation. Then there exists a positive real number XBmwnip) 
such that for almost every Brownian path u starting ati^, j log converges to XBrown(p) 

when t — 7- oo. 

In any holomorphic family (pA)AeA of non-elementary representations of G, we thus obtain a 
Lyapunov exponent function A i— XBrown(pA)5 and we define the bifurcation current associated 
to the Riemann surface structure X as Tbif = dd'\{xBi:ovm{p\))- From now on we will refer 
to it as the natural bifurcation current. We observe that this current depends only (up to a 
multiplicative constant) on the conformal structure on X, and not on the conformal metric 
on X which was chosen to generate the Brownian motion (in the case where X has punctures 
some extra care is needed about the behaviour of the metric near the cusps). Indeed, by 
conformal invariance of the Brownian motion in dimension 2, the Lyapunov exponent function 
associated to any other choice of conformal metric will be a constant multiple of XBrown (see 
Remark 12.81 for details). 

A similar procedure yields a Lyapunov exponent Xgeodesic = 2xBrown by using a generic 
geodesic starting at -k instead of a Brownian path (the factor ^ here is the drift of Brownian 
motion in the hyperbolic plane). In particular Tbif = ^dd'^Xgoodcsic- A statement analogous 
to Proposition |A] in this context was already contained in the work of Bonatti, Gomez-Mont 
and Viana |BGM] (at least when X is compact). 

Yet another way to define the Lyapunov exponent is to use the isomorphism between G 
and a lattice F < PSL(2, M) provided by uniformization. Let IHI be the hyperbolic plane. 
For a sequence r„ increasing sufficiently fast to infinity (e.g. r„ = n° will do), we pick at 
random an element 7„ G Fni?H(0,T„), relative to the counting measure. Then almost surely, 

4^ log l|p(7n)|| converges to XBrown 

As in [DDT] , the first main result asserts that Tbif is indeed a bifurcation current. 

Theorem B. Let X be a hyperbolic Riemann surface of finite type and {px)xeK be a holomor- 
phic family of non- elementary parabolic representations of G = 7ri(A, *). Then the support 
of the natural bifurcation current Tbif is equal to the bifurcation locus. 

The strategy of the proof is to rely on the corresponding support theorem in [DDI]. For this, 
we use a discretization procedure for the Brownian motion, which was devised by Furstenberg 
|F2] and then refined by Lyons and Sullivan |LSj . among others. In some sense, it allows to 
project the Wiener measure (based, say, at *), to a probability measure p. on G. We then 
show that the natural bifurcation current is equal to Tbif^^, up to a multiplicative constant, 
from which the result follows. The point is that the support theorem in j DDlj requires an 
exponential moment condition on p. When X is compact, this fact is apparently part of the 
folklore. We show that when X is merely of finite volume and /> is a parabolic representation, 
p still satisfies the following slightly modified moment property 



I exp (/3(ie(0, /o(g)(0))) dp{g) < oo for some /3 > 0, 
G 
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that is enough to apply the results of |DDlj (see Proposition I2.16p . The proof is rather 
technical, and occupies most of ^ We believe that it could find applications beyond this 
work. 

The second main theme of the paper is the equidistribution of random hypersurfaces of the 
form Z{gn,t) towards the bifurcation current. Given the Riemann surface structure X, there 
are at least two natural choices for the group elements Qn- 

- A first possibility is to pick randomly an element 7^ in F n -Bh(0, r„) (according to 
the counting measure), and let r„ — )• 00. 

- A second choice is to choose at random a closed geodesic 7„ of length at most r„ on 
X. This geodesic defines a conjugacy class in G, so, abusing notation slightly, we can 
define tr^(/5([7])) and it makes sense to consider Z{'yn,t) = {A E A, tr^([7„]) = t}. 

In both cases, we show that the expected equidistribution holds. For concreteness, let us 
state precisely the equidistribution theorem for random closed geodesies (see Theorem 13.11 
below for the other case). The notion of a random closed geodesic will be made precise in 
^3.31 Several natural models are possible, the theorem applies to all of them. 

Theorem C. Let X be a hyperbolic Riemann surface of finite type and {px)x£\ be a holo- 
morphic family of non- elementary parabolic representations of G = 7ri(X, ★) into PSL(2,C). 
Fix t £ C Let (r„) be a sequence such that for every c > 0, the series J2n>o c""^^" converges. 
Let ("fn) be a sequence of independent random closed geodesies on X of length at most rn- 
Then almost surely, 

J— -[Z(7„,t)] Tbif. 

41ength(7„) " n^oo 

Notice that this result does not follow from the equidistribution results in |DD1| since the 
random elements 7„ are not obtained from a random walk on G. To prove the theorem, we 
consider the sequence of psh potentials 

= 71 — V |tr^(PA(7n)) - t\ . 

4 length (7n) 

We need to show that this sequence of psh functions a.s. converges to XBrown(/OA) hi Ll^^{A). 
The main step is to study this convergence for a fixed representation. 

Theorem D. Let X be a hyperbolic Riemann surface of finite type and p : G ^ PSL(2,C) 
be a non- elementary parabolic representation. Let (rn) and (7^) be as in Theorem [0 Then 
almost surely 

n:^^log|tr2(p(7„))| 

XBrown 

41ength(7„) ' ' n->oo 

Notice that when p is the uniformizing representation (i.e. p = id when G is identified to 
F), log |tr^([7„]))| is nothing but the length of 7„ (up to an error of order ©(e-i^'^s*^^^"))). 
This theorem thus describes how this quantity evolves as p is deformed in the space of repre- 
sentations. 

By the Borel-Cantelli lemma, to prove Theorem |D] we need an estimate on the number 
of geodesies of length at most r such that 4iength(7) |fa^(p(7))| - XBrown (/>) > £• This 
counting is based on large deviation estimates for the traces of random products of matrices, 
together with a precise understanding of the Furstenberg-Lyons-Sullivan discretization pro- 
cedure. Observe that the exponential moment property of the discretization measure is also 
crucial here. 
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Finally, Theorem [C] is a rather straightforward consequence of Theorem [Dl by using some 
standard compactness properties of psh functions. 

The outline of the paper is as follows. In after a few facts on hyperbolic geometry and 
PSL(2,C) representations, we recall some necessary results from [DDI]. We also discuss the 
basic properties of Brownian motion on hyperbolic surfaces. In ^ we define the natural 
Lyapunov exponent and prove Theorem |Bj For further reference as well as for completeness, 
we include a detailed treatment of the discretization procedure. The equidistribution theorems 
are established in 33 We conclude the paper with a few open questions. 



The results of this paper were announced in a preliminary form in [DDH §5.1] and Duj 
54.5]. 



1. Preliminaries 

Throughout the paper we use the following notation: if u and v are two real valued func- 
tions, we write u x u (resp u < v) if there exists a constant C > such that < v < Cu 
(resp. u < Cv). We also use the letter C for a generic positive "constant", which may 
change from line to line, independently of a varying parameter which should be clear from 
the context. 



1.1. Plane hyperbolic geometry. To fix notation, we summarize a few basic facts from 
hyperbolic geometry. See Beardon |Beaj for a systematic account. We let (often simply H) 
be the hyperbolic plane (resp. be the hyperbolic 3-space) of curvature —1. The hyperbolic 
distance is denoted by dm- The hyperbolic area in is simply denoted by by Area(-), and 
we recall that the area of a ball of radius r is 47rsinh^ |. As usual, we may switch from 
the Poincare disk (resp. ball) to the upper half plane (resp. space) model depending on the 
particular situation. The group of orientation preserving isometries of (resp. H^) will 
be viewed as PSL(2,M) (resp. PSL(2,C)). The S0(2) (resp. SU{2)) invariant point wih be 
denoted by 0. We let ]]•]] be the operator norm associated to the standard Hermitian metric 
on C^; recall that \\A\\ is the spectral radius of A* A. The norm of an element 7 G PSL(2, C) is 
by definition \\A\\, where A is any matrix representative of 7 in SL(2,C). We also sometimes 

use jj-jjg defined by ]]7]]2 := (Ja)^ + 16]^ + ]c]^ + jd]^)"^'^^, where A = ("^). Likewise, if 
7 € PSL(2,C), tr^(7) and ltr(7)] are well-defined quantities. Norm and displacement are 
related by the identity 

(1) Il7ll2 = 2cosh((ie(0,70)), 

in particular 21og]]7]] ~ dH(0,70) as ]]7]] — 00. If 7 is a hyperbolic (resp. loxodromic) 
element then we denote its translation length by ^(7), which satisfies 

(2) cosh(^) =iltr(7)l. 

1.2. Spaces of representations of 7ri(S') into PSL(2, C). A good source for the material in 
this section is |Kap| . Let be a (topological) surface of finite topological type, that is, homeo- 
morphic to a compact surface of genus g with a finite number n of punctures. Throughout the 
paper we always assume that xi^) = 2 — 2g — n < 0. If S" is endowed with a complex structure, 
the corresponding Riemann surface will be denoted by X. Fix a basepoint -k £ S and put 
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G = 7ri(5, ★). The space of representations of Hom(G, PSL(2, C)) is an affine complex alge- 
braic variety (recall that PSL(2,C) is affine because PSL(2,C) ~ S0(3,C)). A general holo- 
morphic family of representations is simply a holomorphic mapping A — )• Hom(G, PSL(2, C)). 
The group PSL(2,C) acts algebraically by conjugation on Hom(G, PSL(2, C)), and the alge- 
braic quotient Hom(G, PSL(2, C))//PSL(2, C) is the character variety X{G). Coordinates on 
X{G) are given by functions of the form p ^ tr^(/9(7)). We also let X^^{G) C X{G) be the 
subset corresponding to non-elementary representations. 

When S is not compact we actually rather work with the relative character variety Afpar(G), 
which is defined as follows. Let (/ii)i=i...n in G be classes of closed loops turning once around 
the punctures. We define the set of parabolic (or type preserving) representations by 

Hompar(G,PSL(2,C)) = {p G Hom(G, PSL(2, C)), Vi, p{hi) is parabolic} . 

The algebraic quotient ^Vpar(G) = Hompar(G, PSL(2, C))//PSL(2, C) is the relative character 
variety. It can be shown that Afpar(G) is of complex dimension 6(7 — 6 + 2n, and is smooth at 
non-elementary representations. 



1.3. Bifurcation currents. Here we gather some material from [DDlj. See also |Duj 



for 
for 



a general account on bifurcation currents in one-dimensional dynamics, and [Beal Kap 
basics on Kleinian (and more generally Mobius) groups. 

Let G be as above, A be a connected complex manifold and p = (/9a)agA be a holomor- 
phic family of representations of G into PSL(2,C). We usually work under the following 
assumptions: 

(Rl) the family is non-trivial, in the sense that there exists Ai, A2 such that the represen- 
tations p\^, i = 1,2 are not conjugate in PSL(2,C); 
(R2) there exists Aq € A such that p^o is faithful; 
(R3) for every X £ A, px is non-elementary. 

It may sometimes be useful to replace (R3) by the following weaker statement 

(R3') there exists Aq G A such that px^ is non-elementary. 

In that case the subset £^ C A of elementary representations is a proper real-analytic subset. 
According to Sullivan [Su2] there exists a partition A = Bif U Stab of A into an open locus 
Stab of stable representations (in the sense of dynamical systems) and a closed bifurcation 
locus Bif. By definition, a representation p is stable if nearby representations in A are quasi- 
conformally conjugate to p. 

Fix a probability measure p on G, whose support generates G as a semi-group. We will 
assume the following locally uniform exponential moment condition: 

(M) VAo G A, 3iV(Ao) neighbd of Aq, and o" > s.t. VA G iV(Ao), / \\px{g)f dp{g) < 00, 

Jg 

where ||-|| is any matrix norm on PSL(2,C). Under such a moment condition one can define 
a Lyapunov exponent function on A by the formula 

X;.(A)= hm - [ log\\pxig)\\dp^{g), 

and the associated bifurcation current Tbif^^ = dd^x^l ^we deliberately emphasize the depen- 
dence on p here. 

The following result was established in |DDH Thm. 3.13]. 
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Theorem 1.1. Let {G,p,fi) be a holomorphic family of representations ofG, satisfying (Rl, 
R2, R3'), endowed with a measure /i, generating G as a semi- group, and satisfying (M). Then 
Supp(Tbif,^) is equal to the bifurcation locus. 

A comment is in order here: the theorem was stated in [DDI] under the more restrictive 
moment condition 

(3) there exists s > s.t. / exp(s length (5) < cxd, 

Jg 

where the length is relative to some finite system of generators of G. But, as explained in the 
beginning of |DDll §3.3], it still holds under (M). 

It is important to observe that these notions also make sense on (resp. relative) character 
varieties. Indeed let A be a piece of complex submanifold of the set of characters of non 
elementary representations of G into PSL(2,C), A C ;f^^(G) (resp. A C X^J^{G)). Then 
it can be lifted to a holomorphic family of representations, transverse to the foliation by 
conjugacy classes, so we can define a stability (resp. bifurcation) locus on A, which does not 
depend on the choice of lifting. Likewise, the validity of the moment condition (M), as well as 
the value of the Lyapunov exponent, are independent of the lift. Hence there is a well-defined 
bifurcation current on A, and Theorem 11.11 holds . 



1.4. Brownian motion. More details on the material in this paragraph can be found in 
|CC1 IGrii iFL j Let X be a complete hyperbolic Riemann surface, endowed with its hyperbolic 
metric. Brownian motion on X is defined as the Markov process with transition probability 
-Pt(x, •) = pt{x,y)dy, where the volume element dy is relative to the hyperbolic metric, and 
the heat kernel pt satisfies 

d 1 

— pt(x, •) = -Apt{x, •) and lunpt{x, ■)dx = 5x. 

This process is a diffusion with continuous sample paths so it gives rise for every x £ X 
to a probability measure Wx (the Wiener measure) on the space Q-x of continuous paths 
oj : [0, cx)[— )■ X such that w(0) = x (itself endowed with the compact open topology). We also 
denote by the space of continuous paths ci; : [0, oo[— )• X and put W = JWxdx. We often 
the classical notation P^; (resp. Mx) for the probability (resp. expectation) of some event 
relative to Wx- 

It is useful to think of sample paths of the Brownian motion on X ~ r\EI^ as projections of 
Brownian paths on H^. To say it differently, the Wiener measures associated to a point in X 
and of one of its lifts to M.'^ are naturally identified. With our normalization, the linear drift of 
Brownian motion on is equal to 1/2, that is, for Wo-a.e. Brownian path w, dj^{0,uj{t)) ~ | 
as t —7- 00 (see the appendix in |DDp| for an elementary proof). 

We will require the following precise estimate for the radial density k^{t,r) of the law at 
time t of Brownian motion issued from (i.e. Pt{0, •)) in 

(4) k"{t,r) X {l+r-^)-H-\l+r + t/2y\l+r)exp 

(see e.g. [Dal Thm 5.7.2]; notice that the normalization of the heat kernel is different there, 
hence the t/2 instead of t). 



(r - t/2y 
2t 
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2. The Lyapunov exponent associated to Brownian motion 

Let 5 be a topological surface of genus g with n punctures, with a marked point and 
such that x('S') < 0. Recall that G = 7ri(S', In ^2.11 we show that the data of a structure of 
Riemann surface of finite type on S gives rise to a natural Lyapunov exponent function XBrown 
on Hompar(G, PSL(2, C)) (resp. Xpa.riG)). Proving that dd'^XBrown is a genuine bifurcation 
current occupies the remainder of this section. Our main tool is a discretization procedure for 
Brownian motion, which will be discussed thoroughly in ^2.21 From this point, the main issue 
is to obtain estimates for a probabilistic stopping time, associated to the discretization. The 
case where X is compact is substantially simpler, and will be discussed in ^2.31 In ^2.41 we 
deal with the general case. The proof of Theorem 12. 71 is given in ^2.51 In ^2.61 we re- interpret 
these results from the geodesic point of view, from which Proposition [X] follows. 

2.1. Definition. Let X be a structure of Riemann surface of finite type on S (i.e. a neighbor- 
hood of each puncture is biholomorphic to B*). The Euler characteristic assumption implies 
that X is hyperbolic, and we endow X with its hyperbolic metric of curvature —1. Let 
{X,-k) be the universal cover of (X,*). We fix an isomorphism between X and mapping 
* to 0. Then the deck transformation group G becomes naturally isomorphic to a lattice 
r < PSL(2,M). We let 

Pcan bs the (Fuchsian) uniformizing representation of G such that 
r = Pcan{G), and in the sequel we conveniently identify G and F when needed. 

Fix a representation p G Hom(G, PSL(2, C)). If 7 is any loop on X based at and [7] G G 
denotes its class (from now on we drop the square brackets), we may consider the quantity 
log||p(7)ll- 

We are now ready to define the natural Lyapunov exponent of a parabolic representation 
of G, associated to the complex structure on S. Justifying the validity of this definition will 
occupy the remainder of this subsection. 

Definition-Proposition 2.1. Let X be a hyperbolic Riemann surface of finite type and 
p : TTi(X,-k) — )• PSL(2,C) be a parabolic representation 

Define a family of loops as follows: for t > 0, consider a Brownian path uj issued from 
-k, and concatenate w|[q^^] with a shortest geodesic joining uj[t) and -k, thus obtaining a closed 

loop UJf 

Then for a.e. uj the limit 

(5) XBrown(p) = lim \ log ||p {u)t)\\ 

exists and does not depend on uj. This number is by definition the Lyapunov exponent of p, 
associated to the complex structure X . 

Notice that there is a little ambiguity in the definition for there might be several shortest 
paths joining u}{t) and -k. As it will be clear from the proof, this does not affect the existence 
nor the value of the limit in (|^ . It will also follow from the proof that XBrown 

does not depend 

on -k. Since with our normalization, the escape rate of the Brownian motion in is ^, we 
deduce from ^ that XBrown(/Ocan) = 3- 

To prove the existence of XBrown and study its properties, we will need to consider the 
holonomy of non-closed paths on X. For this, we work on the suspension of p. For any 
representation p € IIom(G, PSL(2, C)), the suspension of p is by definition the quotient of 

xP^ under the diagonal action of G. To be specific, Mp = M.'^x¥^ / ~ where the equivalence 
relation is defined by (p, z) ~ {p\ z') if there exists g E G such that {p' , z') = {g{p), p{g){z)). 
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The natural projection vr : Mp — )■ X makes Mp a flat P -bundle over X, and the horizontal 
foliation on X descends to Mp as a holomorphic foliation with leaves transverse to the 
P^ fibers, whose monodromy representation is p. 

The point is to study the growth rate of holonomy between vertical fibers in Mp along 
generic Brownian paths in X, measured with respect to some smooth family of spherical 
metrics on the fibers. When X is compact the choice of the metrics on fibers is harmless. On 
the other hand in the non-compact case this requires some care, so let us dwell a little bit on 
this point. 

By a spherical metric, we mean a conformal metric on P^ of curvature 1. Equivalently, it is 
the Pubini-Study metric associated to some Hermitian structure on C^. Fix a smooth family 
of spherical metrics on the fibers of Mp. If uj is any continuous path on X, we denote by hp{uj) 
its holonomy hp{uj) : 7r~-'-(7(0)) — )• 7r~-'-(7(l)). We define ||/ip(cj)|| as sup2,g^-i(^(o)) IMx/ip(w)||, 
where the norm is relative to the choice of conformal metrics on the fibers. 

We say that the family of spherical metrics is normalized if the induced metric on 7r~^(*) 
is the standard Fubini-Study metric. This means that the natural map {0} x P^ — )■ 7r~^(*) 
induced by the definition of Mp as a quotient, is an isometry when P^ is given the standard 
Pubini-Study metric. This condition is meant to ensure that for closed loops, the definition 
of ||/ip(a;)|| coincides with the natural one. 

By definition a family of spherical metrics on Mp is Lipschitz if there exists a constant 
Cp > such that for every smooth path on X, it holds that 

(6) log||/ip(a;)|| <Cplength(a;), 

where the length is relative to the hyperbolic metric on X. 

When X is compact any smooth family of metrics is of course Lipschitz. In the non compact 
case we have the following result. 

Proposition 2.2. If p is a parabolic representation of G, then there exists a smooth normal- 
ized Lipschitz family of spherical metrics on the fibers of Mp . Furthermore the constant Cp 
in ([6]) can he chosen to be locally uniformly bounded on IIompar(G, PSL(2, C)). 

Notice that the result is false for a general representation. Indeed consider a path of length 
L turning about times around a cusp of X. Let g £ G representing a simple loop aroung 
that cusp. Then if p(g) is loxodromic, we see that ([6]) will be violated. On the other hand it 
can be shown that the proposition (as well as all results in this section) remains true as soon 
as every such p(g) has eigenvalues only on the unit circle. 

Proof. The first observation is that the set of spherical metrics on P^ is parameterized by H^. 
This follows from the fact that such a metric is the visual metric 6{p) from a point p G 
on the sphere at infinity P^. Another way to see it is to observe that a spherical metric is 
defined by a Mobius map from P^ (equipped with the standard Fubini-Study metric) up to 
precomposition by an isometry, and PSL(2, C)/PSU(2) ~ H^. 

We denote by dsp the metric associated to p G H'^. We claim that ii p,q G H^, we have 
that 

(7) log 

This follows from an elementary computation. Assume first that g = 0, so that dSg = ds^ = 
' ^U,2 is the usual spherical metric. Let p £ H^, p = A^^{0) for some A G PSL(2,C), so 



dsp 
ds„ 



dMip,q) + Oil). 
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that dsp = M*dsQ. Putting ^ = we infer that 

\dz\ 



dSr 



\az + + |cz + 



therefore 



dsn 



L4- 



-1||2 



log 



l^ll , and finally 



dsr 



dsr 



21ogp|| =dH(0,p) + O(l). 



To handle the general case, pick an isometry h of bringing q back to 0, and observe that 



since 4^ is a function on P^, 



ds, 



h*dsn 



h*dsn 



r/s 

dSa 



dSr 



dSn 



A smooth family of spherical metrics on a trivial line bundle Z? x is now simply the 
data of a smooth map /:!?—>• H^. A smooth family of spherical metrics on the fibers of 
X P^ descends to a family of spherical metrics on Mp if and only if the associated map / is 
/9-equivariant, i.e. for every g £ G, f o pcanid) = p{g) ° f ■ Furthermore the family is Lipschitz 
iff / is, and the constant Cp is controlled by the Lipschitz constant of /. This follows from 
d?]), since the holonomy of the horizontal foliation on IH^ x P^ is the identity. 

To prove the proposition it is thus enough to show that there exists a smooth Lipschitz 
/O-equivariant map / : H'^ — )• H^. The family of spherical metrics will be normalized provided 
/(O) = 0. To construct such an /, recall that since X is of finite type, there exists a finite 
number of open disjoint subsets Ci, . . . ,Cr oi X such that X \ [_]-Ci is compact, and each 
Ci is isometric to the quotient of {z G | 9z > pi]. by the translation z i— )• z + 1 for 
some positive pi. For every i = l,...,r, we pick some lift Ci C EI^ of Cj, and we first 
define / on the closure of Cj. Denote by G G a generator of the stabilizer of Ci in C. 
By assumption, p{gi) is parabolic, hence conjugate to gi. Thus, there exists an isometric 

immersion f : Ci ^ H'^ such that f o g^ = p{gi) o / (we take the map in PSL(2,C) which 

conjugates gi to p{gi))- We extend the map / by p-equivariance to a map defined on (J^ GCi. 
Then, because is contractible, a classical topological argument shows that the map / can 
be extended as a smooth p-equivariant map from to such that /(O) = 0. We claim 
that this map is Lipschitz. Indeed, is G-invariant (the norm of Df is measured using 

hyperbolic metrics on source and target), thus descends to a function on X. But since / is an 
isometric embedding on Q, we have ||-D/|| = 1 on the cusps Cj's. Hence ||-D/|| is bounded, 
since X \ UiCi is compact, and / is indeed our desired smooth Lipschitz p-equivariant map. 
Finally, it is clear that the Lipschitz constant of / can be chosen to be locally uniformly 
bounded in Hompar(G', PSL(2, C)), so the same is true of the constant Cp. □ 

The following consequence of the previous proposition will be useful. 

Corollary 2.3. For every parabolic representation p there exists a constant fip depending only 
on the constant Cp of ([6]) such that \\p{g)\\ < ||/0can(5)||^''- 
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Proof. Indeed for g £ G, let ag be the geodesic segment [0,pcan(5)0] in H^. Then in Mp, the 
holonomy of the projection of ag in X is p{g)- By ([1]), ^ and the normahzation of the family 
of spherical metrics we infer that 

log \\hp{ag)\\ = log \\p{g)\\ < Cplengthiug) < argcosh ^^^^^^ < Cplog (2 WpcMf) ■ 

Here we use the inequalities II7II2 ^ 2||7||'^ and argcoshx < log 2a;. Hence we get that 
11/0(5)11 < (2 ||/0can(5')||)^'^''- To concludc the proof, just note that since Pcanig) contains no 
elliptic transformations, there exists c > such that ||/9can(5')|| > 1 + c, so there exists a such 
that 2 < II 

Pcan(5') IP- The result follows. D 

We have the following analogue of Definition 12.11 in this context. 

Proposition 2.4. Let X he a hyperbolic Riemann surface of finite type, p : iti{X) — t- 
PSL(2,C) be a parabolic representation, and \\.\\ a smooth normalized Lipschitz family of 
spherical metrics on Mp. With notation as above, for W-a.e. path uj, the limit 

(8) X{p-^^) = lim jlog ||/ip(a;|[o,t])|| 

exists and does not depend on u nor on the the choice of \\-\\. 

Proof. Let a be the shift semi-group acting on by cr<j : uj{-) 1— ?■ u}{- + s). Recall that W = 
J Wxdx. Since the Riemannian measure is invariant under heat flow, W is invariant under 
the time shift Gs : i^(-) 1— s- + s). The family of functionals Ht{io) : to 1— log ||/ip(6<j|[o j]) || 
deflnes a sub-additive cocycle. In order to apply Kingman's subadditive ergodic theorem, we 
need to check that this cocycle is Vl^-integrable, that is, that / log ||(i/ip(a;|[o ^j) || dW{io) < 00 
for some (hence every) t > 0. For this we use the Lipschitz property of the family of spherical 
metrics. Let X be a Dirichlet fundamental domain containing for the action of T on H^. 
For a Brownian path m, let w be its lift to starting from the lift of ct;(0) lying in X. By 
^ we have that 

log||/ip(a;|[o,t])|| < dH(w(0),w(t)) + o!H(0,w(0)). 

Thus 

log 1 1 dhp {uj\[o^t])\\dWiu)< Ex {dm (w (0) , w (t ) ) ) dx 
which is finite thanks to the (super)exponential decay of the heat kernel on hyperbolic plane 

(cf. (BD). □ 



The following result finally justifies the validity of Definition 12.11 

Proposition 2.5. Let X , p, and \\ \\ he as in Proposition \2.4\ Then for W^, a.e. Brownian 
path u), we have that 

lim jlog||p(a;t)|| = lim j log ||/ip(a;|fo,ti) || = XBrown(p)- 

Furthermore the limit in ([8]) is also equal to XBrown(p)- 

Proof. Let us first show that for VFj, a.e. uj, lim(_i,oo 7 log ||/ip(u;|[o^i]) || exists and has the 
same value as the one in ([8]). Indeed the diffusion of the Dirac mass (5^ at time 1 is abso- 
lutely continuous with respect to volume measure on X. Now we simply write /ip(u;|[o^t]) = 
^p('^l[i,t]) ° ^p('^l[o,i]) aiid by the Markov property of Brownian motion together with Propo- 
sition [231 we get that for VK* a.e. oj, 7 log ||w|[i^t]|| converges to the limit in ([8]). 
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To finish the proof, we show that for a.e. cj, 



(9) 



log \\p{u;t)\\ - log ||/ip(a;|[o,t]) 




0. 



If X is compact this is obvious since dx{u}{t),-k) is uniformly bounded. In the general case, by 
the Lipschitz property we have to show that for W^, a.e. u, dxiu}{t),-k) = o{t). A way to see 
this is to use the following two facts: first, a theorem of Sullivan [Sulj asserts that for almost 
every v in the unit tangent space Si,X to the unit speed geodesic ray 7(t) issued from 
-k pointing in the direction v satisfies (ix(*,7(i)) = 0{\ogt). The other result we use is the 
classical fact that Brownian paths on hyperbolic space are shadowed by geodesic rays (see e.g. 
I An] ) . More precisely, if w is a Brownian path starting from G H^, then limt-^oo^(i) = ^oo 
exists in d'M?, and if ^x,u)aa denotes the unit speed geodesic ray joining x and Wqo, then for 
Wo a.e. cj, d^{ijj{t),^x,u)oa{\'^) — o{t) (this is actually O(logt) by a theorem of Ancona [AnJ). 
Note that the distribution of Wqo is absolutely continuous on the circle at infinity (with density 
given by the Poisson kernel). Therefore on X, for a.e. a;, there exists a unit tangent vector 
such that dx (<^(^)) 7*,t; (^*)) = o{t). Since the distribution of the tangent vector v is induced 
by that of it is absolutely continuous with respect to Lebesgue measure on S^X. So we 
can combine this estimate with Sullivan's, and conclude that (ix('^(i),*) = o{t) a.s., which 
imphes ©. □ 

For further reference, we now make the easy observation that the natural Lyapunov expo- 
nent is insensitive to finite coverings. 

Proposition 2.6. Let X he a hyperbolic Riemann surface of finite type, and vr : (Y^yo) — )• 
(X, he a finite covering. Any parabolic representation p : 7ri(X, ★) — )• PSL(2,C) induces by 
restriction a parabolic representation py of niiY^yo). Then, denoting by XBrown 

{PY,Y) the 

Lyapunov exponent of py associated to the Riemann surface structure on Y , we have that 

XBrown — XBrown 

Proof. This is just a matter of following step by step the proof of the existence of XBrown(P5 X)- 
Indeed, we observe that the suspension Yp^ is naturally a covering of Mp, and the lift of a 
Lipschitz metric ||-|| on the fibers of Mp is a Lipschitz metric on the fibers of Yp^. As already 
observed, the Wiener spaces W-^ on X and Wy^ are naturally identified to Wq in 'M? . Now 
by construction, if w is a Brownian path issued from t/q in Y, we have that ||/ipy (a;|[o^t])|| = 
||/ip(7r(a;)|[o,<])||- The result follows. □ 

Let now {p\)\^\ be a holomorphic family of parabolic representations of G. Thanks to 
Definition 12.11 we can define a Lyapunov exponent function XBrown on A, and we put Tbif = 
XBrown- We will refer to it as the natural bifurcation current of the family. The following 
result asserts that Tbif is indeed a bifurcation current. 

Theorem 2.7. Let X be a hyperbolic Riemann surface of finite type, and {p\)x^\ be a holo- 
morphic family of representations of G = iti{X), satisfying (Rl, R2, R3'). Then XBrown is 
a nonnegative psh function on A, and the positive closed current Tbif = '^c^'^XBrown satisfies 



If furthermore (R3) holds then XBrown is positive and locally Holder continuous. 

The proof consists in exhibiting a measure ponG satisfying (M) and such that Tbif^p = Tbif 
(up to some multiplicative constant). It then suffices to apply Theorem ll.il It will be carried 
out in the next subsections. 



Supp(Tbif) = Bif . 
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Remark 2.8. We have constructed the natural bifurcation current by using the hyperbohc 
metric on X. Let us observe that when X is compact any other conformal metric would 
essentially yield the same object. Indeed, choosing a different conformal metric on X only 
changes the "clock" of Brownian motion, that is, sample Brownian paths for the new metric 
are time-reparameterized Brownian paths relative to the original one. It then follows from 
an easy ergodicity argument that the associated Lyapunov exponent function is a constant 
multiple of XBrown- When X has punctures, the same statement holds, provided one restricts 
to an appropriate class of conformal metrics, for instance the ones that coincide with the 
Poincare metric outside some compact set of X. 

2.2. Discretization of the Brownian motion. We will rely on a discretization procedure 
for the Brownian motion which was introduced by Furstenberg |F2] . and subsequently stud- 
ied by several authors, in particular Lyons and Sullivan [LSj (see also Kaimanovich [Kaij . 
Ballman-Ledrappier |BaLej . Ancona |Anj ). Since we will need precise moment estimates, we 
include a detailed treatment of the discretization, which is close (but not identical to) that of 
Lyons and Sullivan. 

Let r be a finite covolume lattice in PSL(2,IR), and X = F \ M. Let G M and ii > 
such that the balls 3(^0, R) are disjoint, for 7 G F. Since the Brownian motion on X is 
recurrent |Gri| . for < r < i? the set dr = \J^dB('yO,r) is recurrent, namely a.e. Brownian 
path starting at a point x hits dj. in finite time. We denote by (resp. Tr) the first hitting 
time to dr (resp. dji). Since dr (resp. dji) is closed, Tr (resp. Tr) is a stopping time. 

2.2.1. Sub-balayage. Recall that the balayage of a finite measure m on a closed recurrent set 
E is defined to be the image of the measure J Wxdm{x) under the map oj i— )• oj{Te) where 
T = Te{u}) is the first moment where uj hits E. 

We denote by Bp, Sp respectively the ball and sphere centered at and of radius p. If y is 
a point interior to Bji, the balayage of the Dirac mass at the point x on is a probability 
measure absolutely continuous with respect to the spherical measure ds on Sr (conveniently 
normalized so that ds is a probability measure), the density being given by the Poisson kernel. 
Hence, there is a constant < p < 1, such that if x is chosen in the sphere Sr, its balayage 
on Sr is bounded from below by pds. 

The following lemma will allow us to select a "good" set of paths from Sr to Sr. For every 
X G *S'^, a. e. Brownian path issued from x intersects Sr. The hitting time is Tr(uj) and we 
denote by e(a;) = uj(Tr{uj)) the exit point. We use the standard notation {J-'t)t>o for the 
Brownian filtration. 

Lemma 2.9. For every x G Sr there exists a Txji-measurahle set Gx, with P(Ga;) = p such 
that e^,{lG,^Wx) = pds. 

In addition, reducing p slightly if necessary, Gx can he chosen so that there exists a constant 
M such that for every x £ Sr and lo G Gx, Tr{u;) < M. 

Proof. Let us first establish the lemma without the bound M on the exit time. Denote by 
bx{y)ds{y) the balayage measure of 6x on Sr. Notice that bx{y) > p for every y G Sr. For 
< t < 1 consider the piece of circle 



We use Ct as a "continuously opening barrier" for Brownian paths going from x to Sr. For 
< t < 1, let £^t C be the set of paths starting at x that avoid Ct before hitting Sr. For 
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t = the barrier is closed so = while for t = 1 it is open and Wx{Et) = 1; observe also 
that for t < t' , Ct Z) Cf so Et C Ef. Therefore, for every y G Sr, fy'.t^ Wx{Et\e~^ {y)) is 
a continuous increasing function with /j/(0) = and fy(X) = 1- The continuity of fy follows 
from the fact that the set of Brownian paths (starting from x and exiting at y), intersecting 
Efi \ Et has small measure when t' — t\s small. Let t{y) G (0, 1) be the unique value of t such 
that fy{t{y)) = 5^^- We can now define 

Gx = {uj e Clx, i^\[0,Tii{oj)] n Ct{e{u})) = 0} ■ 

Clearly, this set is Tt^ measurable (recall that e(a;) = u!{T}i{uj))). Since by definition Gx fl 
e~^(y) = ^t{y) l~l ^~^{y) we infer that 

e.{lG.Wx) = J Wx{Gx\e-\e{u)))dWx{uj) = J Wx{Et^c(u.)))\e~\e{uj)))dWxito) 

= WxiEt(^y-j\e~^{y)) e^Wx = Wx{Et^y)\e~^{y))bx{y)ds{y) = pds, 

which was the desired property. 

To prove the second claim, observe that for every pi < p, if Ax{M) C Wx denotes the set 
of paths such that Tr(w) < M, then for M sufficiently large, e^,{lA^(Tl)Wx) > pids, uniformly 
in X -this is a consequence of Dini's Theorem. Prom this point, to get the result it suffices to 
repeat the above proof with l^^^^^^Wj, instead of Wx- □ 

We extend by equivariance the definition of the sets Gx to all of dr by setting G^x = jGx ■ 

2.2.2. A stopping time. For every point xq G H and every lo G ^xq^ '^^ define a sequence of 
stopping times 

Uo = Tniuj) <Ui= Ui{u}) < Ui = Ui{uj) < . . . < Uniu}) < Uniuj) < ■ 

as follows: Uq = Tji{uj), and by induction, Un is the first time after Un-i at which u hits dr, 
that is Un{(^) = Un-i{uj) + Tj.{aur,-i{'^)), where as before at denotes time shift by t. Likewise, 
Un is the first time after n„ at which co hits Or, i.e. Un{uj) = Un{oj) + Tii{au„{^))- 

By the strong Markov property of Brownian motion, almost surely there exists an integer 
n such that ctu^w belongs to G^^^^y We define T{uj) := Uk{uj), where k = k{oj) is the first 
such integer. The measurability property of the sets Gx ensures that T is a stopping time 
(intuitively at time t we know if we already went from dr to Or through some Gx)- We let 
xi{(jj) = (jj{T{u)) and ri(u;) be the element of V such that x\ G dB{ri ■ 0,R)- The law of xi 
defines a transition probability n(a;o, •)• construction it satisfies the equivariance property 
n(7Xo, •) = 7*n(xo, •) for 7 G r. 

Repeating this process independently gives rise to sequences {xn)n>o € (H^)^ and (r„)„>i G 

2.2.3. The discretization measure- It follows from the construction that whatever the distri- 
bution of the point xq is, the distribution of the point xi is a convex combination of the 
spherical measures ds^ = 7*ds on Or- Hence, if xq has distribution ds on Sr (observe that 
C/q = in this case), then xi has distribution ds for a certain measure on T. 

We now define a Markov chain (x„) with initial distribution given by the normalized spher- 
ical measure ds on Sr and transition kernel 11. By the equivariance property of 11, the 
distribution of x„ is f/^ * ds- 
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Furthermore, since the transition kernel is defined by stopping the Brownian motion, by the 
strong Markov property we infer that the map w i— )• (r„(a;)) sends the measure J^^ Wxds{x) 

to the path measure P on of the right random walk on F induced by /j. 
By definition fi is the discretization measure (notice that it is not unique). 

2.3. The cocompact case. The core of the proof of Theorem 12.71 is the fact that the dis- 
cretization measure satisfies the exponential moment condition (M). When X is compact, this 
fact is well-known to the specialists (for instance it is announced without proof in [Anj ) . but 
also not so easy to find in print. Let us explain this case first. Recall that F = /Ocan(G'). 

Proposition 2.10. // the lattice F is cocompact, then the discretization measure jj, satisfies 
the exponential moment condition ^ with respect to word length in G. 

Proof. By the stronger assertion of Lemma 12.91 for every pi < p there exists a deterministic 
time Ti such that for every x £ Sr there exists a measurable set Gx such that e^{lc^Wx) = 
pids and moreover for every oj £ G]., Tr{oj) < ti. 

On the other hand, since X = F \ H is compact, there exists < p2 < 1 and T2 > such 
that for every x E H, the probability that a Brownian path starting at x hits dr in time less 
than T2 is bounded from below by p2- 

Therefore, by the Markov property we infer that for every x E EI and every to, 

Px(T > to + Tl + T2) < (1 - P1P2) ■ fx{T > t), 

where T is the stopping time defined in ^2.2.21 From this we deduce that there exist constants 
C,a > independent of x such that 

(10) ¥x{T >t)< Ce-"*. 

Since the discretization measure is the image of J^^ Wxds{x) under to 1— )• uj{T{uj)), and the 
distance between ri{uj) ■ and uj{T{lo)) is R, for every 6 > 0, we have that 

/i(7 E F, d{jO,0) > (5) < y ¥xid{u}{T{u})),0) > 5 - R) ds{x), 

where d{-,-) denotes hyperbolic distance. Fix > 1 as in Lemma 12.111 below. For every 
X £ Sr, we have that 

¥x{d{uj{T),0) >6-R)< Fx{d{uj{T),x) > 5 - 2R) 

<¥x\ sup d(a;(t),x) > (5-2/?) ) (r> ^-^^ 

\t<(S~2K)/K J \ K 

By Lemma [2.11l below the first term decreases exponentially fast when 6 tends to infinity, while 
the second does because of (fTO]l . By the homogeneity of hyperbolic plane, these estimates are 
uniform for x E Sr, and me conclude that /i(7 E F, ci(7(0),0) > 6) decreases exponentially 
with S. Finally, since F = yOcan(G) is cocompact, (FO,d) is quasi-isometric to any Cayley 
graph of G and we are done. □ 

The following result is presumably well-known but we could not locate a precise reference. 

Lemma 2.11. There exists a positive constant K such that 



sup d(0,a;(s)) > Kt 

0<s<t 



decreases exponentially with t. 



16 



BERTRAND DEROIN AND ROMAIN DUJARDIN 



Proof. Let us work in the disk model of hyperbolic plane. The expression for the hyperbolic 
distance is exp((i(0,z)) = j^j^j . Since x i— )• is convex and increasing, z i— )• exp((i(0, z)) is 
a subharmonic function. Prom this we infer that if uj{t) is the Brownian motion starting at 
0, the process defined by ( exp((i(0, is a positive submartingale. 

From the asymptotic estimate for the heat kernel given in gH we infer that Eq ( exp{d{0, uj{t)))) 
is finite for every t. Now we have that 

(11) Eo(exp(d(0,a;(t)))) < Eo(exp(d(0, w([t])))) by the submartingale property 

< Eq exp d{uj{i),ui{i + 1)) I I by the triangle inequality 

\ \^=^ II 

= Eo( exp((i(0, a;(l)))) by the Markov property 
<Eo(exp(d(0,w(l))))*+' 

(for the equality on the third line we take successive conditional expectations and use the ho- 
mogeneity of the hyperbolic plane). The Doob martingale inequality applied to exp(d(0, ui{t))) 
asserts that for every a > 0, 

(12) Po f sup expd(0,w(s)) > a ) < -Eo(exp(d(0,tj(t)))) 

\o<s<t ) a ^ 

Hence if K is any real number larger than Eo(exp(d(0,cj(l)))) , applying inequality (jl2p with 
a = ey.'p{Kt) together with ()lip finishes the proof. □ 

Remark 2.12. It is likely that the result holds for any K > ^. K possible way to achieve this 
would be to refine the subdivision of the interval [0,t] in (|lip . 

2.4. The finite volume case. We will need the following estimate on the hitting time of 
the Brownian motion on a given subset of X. We thank J. Franchi for explaining it to us. 

Proposition 2.13. Let X he a hyperbolic Riemann surface of finite type, endowed with its 
Poincare metric, and Liouville measure dx. Let K C X be a compact set of positive area and 
define a stopping time for the Brownian motion by 

Tk{oj) = inf {t > 0, oj{t) G K} . 

Then there exists a constant a > such that for every x G X, 'Ex{e°''^'^) < +oo. Furthermore, 
the function x i— ?■ Ea;(e"^^) < +oo is locally bounded and in L'^{X,dx). 

Of course this uniformity is highly non-trivial for some paths go deeply inside a cusp before 
reaching K. 

Proof. Let (n()t>o be the Brownian semigroup on X, associated to ^A. The spectral theory of 
the Laplacian on X implies that Ilf is strongly contracting, specifically there exists a constant 
c > such that for every / G L'^{X) of zero mean, ||nf/||^2 < e"'^* 11/11^2 (see e.g. [FLJ for 
an elementary approach to the Poincare inequality on X, which classically implies a positive 
lower bound for the first eigenvalue of the Laplacian). Therefore the family of operators 
defined by Qt = y /q converges strongly to / i— )• fxf = f f{x)dx in L^. More precisely 
there exists a function c(T) with lim+oo c(T) = such that 

Qrf- [ f <c{T)\\f\\^,. 

JX L2 
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Indeed write / = /o + // with / /o = 0. Notice that ||/o||i2 < 2 ||/||j^2. Then we have that 



Qrf- f 



L2 



WQtM 



L2 



1 1 

To 



2 

< — 
- ^ 



L2 



+ 



T-T( 



L2 



In this setting, a theorem by Carmona and Klein jCKi Thm 1] asserts that x i— ?■ 'Ex{e°''^^) 
belongs to L'^{X), in particular it is finite a.e. Furthermore, elaborating on [CKl Rmk 1], 
we see that Kx{e°''^'^) is locally bounded. Indeed let to > 0, and use the Markov property to 
write 



J{n<t} Jin>tn\ J 



'{Tk>to} 



Il2 , 



where Ilt{x,y)dy is the transition kernel. Now, as it easily follows from the estimates in 3131 
for fixed to the heat kernel ptoix, •) is in L^, locally uniformly with respect to x. The result 
follows. □ 



From this proposition we get estimates on the stopping times Tr and T/j in H. 

Corollary 2.14. There exists uniform constants aQ and M such that for every x G dn, 
Ej;(exp ao^^r) < (resp. for every x G dr, Ea;(expaoTR) < M). 

Proof. To prove the corollary we work on X. Abusing notation we set dr = nSr, and similarly 
for Or, where vr : — )■ X is the natural projection. Recall that by assumption B{0,R) 
projects to a ball on X. Assume that x £ Or, and fatten dr on the side opposite to x to make 
it a closed annulus dr- Then the hitting times of dr and dr are the same, and we conclude by 
Proposition 12.131 The other case is similar. □ 

The next step is to show that the stopping time T of l2.2.2l defining the discretization admits 
an exponential moment. 

Proposition 2.15. With notation as in there exist constants a and M such that for 

every x £ Sr, 

E^(exp (aT)) < M. 

Proof. Let x G Sr. Recall that for a.e. u G Wx we defined in ^2.2.21 the integer k(u}) as the 
first time at which u}{un{uj)) enters the good set G^(^u„)- For a constant a to be determined 
later, we decompose the integral / eyi'p[aT)dWx as a sum 

V / exp(ar)dW, 

where Vl^ = {w | k{(jj) = K}. We define a sequence of shifted paths w;, < / < 2K — 1 by 
Wo = w, wi = au^{ijj), i02 = cr;7^((jj), etc. We also define for 1 < / < 2K — 1, yiioj) = u;i(0), 
so that for even I (resp. odd /), yi € Sr (resp. yi G Sr). Notice that by definition of K, for 
l = 2j - 1, j < K,uji^ while uj2K-i G Gy^j^_-^. 

Writing T{u}) = Tr{ujQ) + Tr(u;i) + • • • + Tr{u}2K-i), taking successive conditional expec- 
tations with respect to the events T,(uji) = (with • = r or R) and applying the strong 
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Markov property yields 

(13) [ exp{aT)dW^= [ exp{aTriojo))dW^{uJo) ■ ■ ■ 

••• / exp{aTii{uj2j-i))dWy2j_-,{u}2j~i) / exp{aTr{u)2j))dWy2j{uj2j) 

••• / exp{aTR{u}2K-i))dWy^^_-^{u2K~i)- 



Now we claim that for every /3 > there exists a > such that putting q = 1 — p we have 
(i) yx£ Sr, J ex.p{aTr)dW^ < 



(ii) Vy G Sr, Jqc ex.p{aTfi)dWy < qe^ and ex.p{aTji)dWy < pe^ . 
From this and ()13p we immediately get that 

(14) / exp{aT)dW^ < e^^^q^'^p. 

Thus, choosing /? such that e^^q < 1 and summing over K, we obtain that 

/ exp{aT)dW^ < ^ ^(e^^^g)^ < +oo, 
J 9 ^ 

and the result follows. 

It remains to prove our claim. We only deal with (jn]), the other case being analogous. Let 
< a < ao/2, where ao is in Corollary 12.141 let = (3/2 and set E = Gy or E = G^. 
Applying the Chebychev inequality, Corollary 12.141 implies that for every y £ dr, we have 

Wy{TR >t)< Mexp(-aoi). 

We partition the set E as 



E= {En {aTn < /3'}) U [J {E n + an < Tr < ^' + a{n + 1)}) 

n=0 

and decompose the integral exp{aTR)dWy accordingly, which yields 
j^exp{aTR)dWy < e^Wy{E) + - ^ + ") 

<e/^V.(i?) + Cexp (-^) 

with C = ^ ° -aovi • Note that Wx{E) equals p or q, that are both fixed positive numbers, so 
that if Q is sufficiently small we have that 

j^exp{aTR)dWy < e^'W^iE) + Cexp < e^W,{E). 

and claim dn]) follows. The proof is complete. □ 
We can now show that the discretization measure satisfies the moment condition (M). 
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Proposition 2.16. Let X be a hyperbolic Riemann surface of finite type, p be a parabolic 
representation of G = tti{X) and fj, be the discretization of the Brownian motion. Then there 
is a constant a = ap > such that 

(15) I \\p{g)rdp{g) <oo. 

Furthermore, the constant Op can be chosen to be locally uniform in Hompar(r, PSL(2, C)). 
In this case the value of the integral in (|15p is also locally uniform. 

Proof. For p = Pca,n, given Proposition 12 . 1 5l the proof is identical to that of Proposition 12.101 
However in this case we cannot transfer the moment condition in F to a moment condition 
with respect to word length in G (see Remark 12.171 below). Instead we use the Lipschitz 
property of the spherical metrics on the fibers of Mp to obtain moment estimates for p{G). 

Let us be more precise. The proof of Proposition I2.1UI shows that the exponential decay 
property pil|) implies that 

M7er, d(7(0),0) >5) <e-^^ 

for some positive constant c. Since ||7|| x logd(70, 0) we infer that there exists a > such 
that 

{g)\\'^ dp(g) < oo. 

Applying Corollarv 12.31 thus finishes the proof. □ 

Remark 2.17. When X is not compact, the moment condition in (|15p cannot be replaced by 
([3]). Actually, p does not even have a finite moment with respect to word length in G, i.e. 

length(g)dp{g) = oo. This was shown (somehow implicitly) by Guivarc'h and Le Jan in 
[GL] . See also }DKNl Corollary 1.22] for a different approach, as well as Gruet [Gru, Rmk p. 
500] for a neat computation in the case of the modular surface. 

2.5. Proof of Theorem 12.71 Let p be the discretization measure and x^l be the associated 
Lyapunov exponent. We showed in Proposition 12.16) that p satisfies the moment condition 
(M). Furthermore, it is easy to see that Supp(/i) = G. Indeed, recall the balayage construction 
of p. We start with the circle measure ds on Sr and "sub-balayage" it successively on 
Or = T ■ Sr. The first step of the process is to balayage ds onto dj- = T ■ Sr- It is clear that 
the image measure has positive mass on every circle 75*^, 7 S F. Now for each circle jSr, 
there is a set of positive measure of good paths, that leaves some measure on 'jSr (for the 
remaining paths, we repeat the process over again). Therefore p puts some mass on {7}. 

By Theorem 11.11 Supp{dd^Xfi) = Bif. To prove the theorem, we show that there exists a 
positive constant c such that XBrown = cxp, hence Supp(Tbif) = Supp(Tbif,^) = Bif. Under 
(R3), the positivity and Holder continuity properties of XBrown then follow from classical 
results due to Furstenberg |Flj and Le Page |Le) . 

Let T be the stopping time defining p, and T„ be the sequence of times obtained by 
repeating the stopping process. Let Wsj^ be defined by Wsj^ = Jg^Wxds{x). For Wsj^ a.e. 
{Tn) G rn{oj)SR, with (r„) G F^ as in E2l 

Identify F and G through /Jcan and consider a general parabolic representation p. Let as 
before (Tr„ denote the geodesic segment joining and r„(0). Since ci(r„(0), a;(r„)) = i? is 
constant we infer that 

- |log ||/ip(^.„H)|| - log \\hp{uj{Tn))\\\ = - |log \\p{rn)\\ - log \\hp{uj{TM 0. 



r 
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Since the law of r„ is /i" we thus deduce that 

(16) = lim i / log\\p{g)\\dfi^{g) = lim - [ log ||/ip(u;(r„))|| dT^s«(a;). 

n-5>oo n Jq rt -s-oo n J 

To conclude the proof, it remains to show that there exists a constant r > such that — — >■ r 
for a.e. uj. This, together with Proposition 12.51 implies that the right hand side of ()16p 
equals rxsrown 

To prove this, we view Sr as a subset of X. Let 6't be the time shift by T{uj), acting on the 
set of Brownian paths uj issued from Sr, equipped with the measure Ws^- By construction the 
law of ui{T) is ds{x). Prom this fact and the strong Markov property we see that {9t)*Wsii = 
Wsn- Furthermore there cannot be any measurable ^^-invariant subset, again by the strong 
Markov property together with the fact that for every x G Sr, the law of uj{T), conditioned 
to start at x, is ds{x). Therefore 9t is W^^-ergodic 

To conclude, observe that Tn{uj) = J^^Zo'^i^T^) ^^'^ recall from Proposition 12.151 that 
T is integrable. Thus, property (ii) follows from Birkhoff's ergodic theorem. The proof is 
complete. □ 

2.6. The Lyapunov exponent associated to geodesic flow. As before let X be a Rie- 
mann surface of finite type with a marked point endowed with its hyperbolic metric. Por 
V G TJ^X a unit tangent vector, we let 7„ be the associated unit speed geodesic ray. These 
geodesies form another natural class of paths on X. We also denote by ^v,t the loop obtained 
by appending to 7^|[o,t] a shortest path returning to * (as before, we disregard the possible am- 
biguities happening here). This point of view allows to define a Lyapunov exponent function, 
similarly to Definition 12.11 It turns out that it coincides with XBrown up to a multiplicative 
constant. Proposition lAl immediately follows. 

Theorem 2.18. Let X be a hyperbolic Riemann surface of finite type, endowed with its 
Poincare metric. Let p : 7ri(X, *) — PSL(2,C) be a parabolic representation, and \\.\\ a 
smooth Lip schitz family of spherical metrics on the suspension Mp. Then for a.e. v the limit 

Xgeodesic(/?) := ^lim ^ log \\hp{-fv\[o^t]) II = ^lim ^ log \\hp{%t) \\ 
exists and is equal to 2xBrown(p) (in particular it does not depend on v). 

We leave the reader check that the same result holds for geodesies rays issued from generic 
{x,v) G T^X (relative to the Liouville measure), that is, lim(„i.oo 7 log ||/ip(7(2:,t,)|[o,t]) || = 

rown 

Proof. The work was already done in Proposition 12.51 It follows from ([6]) that if uji and 102 
are two paths on X with uji{0) = uj2{0) and such that the respective lifts uji and uj2 to H^, 
starting from the same point, satisfy dMiuji{t),uj2{t)) = o{t), then 

lim y |log||/ip(wi(t))|| -log||/ip(L^2(t))||| =0. 

As already said, generic Brownian paths on are shadowed by geodesic rays. That is, if 
w is a Brownian path starting from G H^, then limt_j.oo w(t) = oJoo exists in dM?, and 
de(a;(t), 7o,a;oo (i^)) = c»(t), where if 7o,tjoo Is the unit speed geodesic ray joining and Wqo- 

To prove the theorem, we use this reasoning in the reverse direction. Taking endpoints 
of semi- infinite geodesies gives the natural "visibility" identification between TqX and 
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sending circle measure on T^X to the harmonic measure associated to on dM."^. For 
V £ Tlx, let uJoo{v) be the associated endpoint. For a.e. Brownian path uj conditioned to 
converge to uJoo{v) we have that 

hm - \\og\\hp{^,{t/2))\\-\og\\hp{u{t))\\\=Q. 

t—^oo t 

Now by Proposition 12.41 for W^, a.e. uj we have that lim j log ||/ip(a;(t))|| = XBrown- Hence 
lim J log ||/ip(7^(t))|| = 2xBrown(/3)- To deal with the loops ^y^^t, as in Proposition 12.51 we use 
Sullivan's result that almost surely, dx(7t)(t),*) = 0{logt). The proof is complete. □ 

Remark 2.19. In |BGM] . Bonatti, Gomez-Mont and Viana study PSL(n, C) cocycles over 
hyperbolic and symbolic dynamical systems, with applications to suspensions of fundamental 
groups of compact surfaces. In particular they show, that if p : Tri{S) — )• PSL(2,C) is a non- 
elementary representation, then the foliated geodesic flow of Mp admits non-zero Lyapunov 
exponents (i.e. Xgcodcsic(p) > 0). We see that the above methods (discretization of Brownian 
motion, etc.) give a new approach to this result, as it was recently observed by Alvarez [Al] 
(for compact S). The novelty here is that we are able to deal with parabolic representations 
on Riemann surfaces of finite type as well. 

3. EQUIDISTRIBUTION in PARAMETER SPACE 

Let as before X be a structure of Riemann surface of finite type on S. Identifying {X,-k) 
with (H, 0) makes G = 7ri(5, isomorphic to a lattice F < PSL(2,M). From now on, we 
consistently identify G and F. We endow F with the metric structure induced by EI on F • 0, 
that is for 7 G F we put ^(7) = de(0,70); observe that ^(7) ~ 2 log ||7||. Notice that when 
X is not compact, this distance is not quasi-isometric to word-length in G. We also put 
Br{r) = {jGT, ^(7) < r)}. 

This metric structure gives rise to a natural notion of a random diverging sequence in F: fix 
a sequence {rn)n>i increasing to infinity and independently choose 7^ G Br{rn) relative to the 
counting measure. For technical reasons we will assume that the series Yln>o converges 
for every c > 0. The resulting random sequences (7^) will be qualified as admissible. 

Let now A be a holomorphic family of parabolic representations of G into PSL(2, C) (resp. 
a holomorphic family of representations modulo conjugacy) satisfying (Rl, R2, R3). For 
t G C, we consider the subvariety in A defined by 

Z(7,t) = {AGA, tr\p^{^)) = t} 

(which is viewed as a variety, that is, possibly with multiplicity), and the associated integration 
current [Z(7,t)]. By convention if Z{'y,t) = A, we put [Z(7,f)] = 0. 
Our purpose in this section is to establish the following theorem. 

Theorem 3.1. Let T be a torsion-free lattice in PSL(2,M), and {px)x£A be a holomorphic 
family of parabolic representations of T into PSL(2,C) satisfying (Rl, R2, R3). Let Teif = 
dd" 

XBrown ^^^6 natural bifurcation current, associated to Brownian motion on X — F \ H. 
// (7^) is an admissible random sequence in T , then almost surely the following convergence 
holds: 

1 

— — — [Z(7„,t)] — > Tbif. 

4d(7n ) ri^oo 
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As already said in the introduction, this result cannot be deduced from the equidistri- 
bution results in [DDlJ. The main step of the proof, achieved in ^3.11 is to show that 
if p G Hompar(r, PSL(2, C)) is a fixed representation, then the exponential growth rate of 
tr^(/9(7„)) is a.s. given by XBrown(p)- This is based on the detailed understanding of the 
asymptotic properties of random product of matrices associated to discretized Brownian mo- 
tion. 

Then in ^3.21 we use Fubini's theorem and subharmonicity to get the almost sure L^^^ 
convergence of A i— >• log |tr^(/>;^(7„)) — t|. Prom this. Theorem 13.11 immediately follows. This 
simple argument was suggested to us by D. Chafai. Notice that the same approach leads to 
a simplification of the almost sure equidistribution theorems in [DDlj as well (see Theorems 
B and 4.1 there). 

In ^3.31 we show that this result can be nicely interpreted in terms of the holonomies of 
random sequences of closed geodesies on X (Theorems |D] and . 



3.1. Distribution of traces in p{T). Let T be a lattice in PSL(2,M). If T is torsion-free 
and yO is a non-elementary type-preserving representation of P into PSL(2,C), the results of 
^enable to define the natural Lyapunov exponent XBrown(/o)- If now P has torsion, by the 
Selberg Lemma, we pick a finite index torsion-free subgroup P' < P, and define XBrown(p) to 
be equal to XBrown(/5|r')- Proposition 12.61 the value of the Lyapunov exponent does not 
depend on the choice of P'. 

Our main result in this paragraph is the following. 

Theorem 3.2. Let T be a lattice in PSL(2,M), and p be a non- elementary type preserving 
representation of P. Then for every e > there exists c{e) > such that for large enough 
r > 0, 



# 7Gi?r(r), 



> e ^ < e 



(l-c(e))r 



^log|tr^(/j(7))| -XBrown(p) 

The delicate point in this theorem is the following: several times so far we have used the 
idea that closing a long geodesic or Brownian path by some path of bounded or short length 
does not affect much the norm of the holonomy. There is no such result for the trace: if 
7 (resp. a) is a "large" (resp. "small") element in P, tr^(p(7)) may be very different from 
tr^(/j(7a)) (even for the identity representation). 

The idea of the proof is to rely on the very precise understanding of generic random products 
7i • • • 7n) where the 7j are chosen according to the discretization measure /i defined in ^2.2[ 
In Proposition 13.31 we show that we can perturb these elements in p{T) by keeping track of 
the traces. We then use the properties of Brownian motion to compare and the counting 
measure, so as to cover most of a large ball in P. 



Throughout the proof c{e) stands for a quantity which may change from line to line, 
dependent on e (and p) but not on n. 

Our first result does not depend on Brownian motion. 

Proposition 3.3. Let P be as above and p : P — )• PSL(2, C) be a non- elementary represen- 
tation. Fix a probability measure p on T, generating T as a semi-group, and satisfying the 
exponential moment condition Jp \\p{'y)\\°' dp{'y) < oo for some a > 0. 
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Then there exists a constant c 
enough 



(17) 



7, 3a G r s.t. \\a\\ < e 



c{n, p) > such that for every e > 0, when n is 
1 



and 



n 



log|trp(7a)| - x^l 



> e > < exp(— cen) 



We note that in this statement we could replace F by an abstract finitely generated group, 
by replacing the condition "||a|| < e*^^"" in (|17|) by "length(a) < cen". 

Proof. Without the perturbative term a, this was proven in |DDH Cor. A. 2] (see also [ Aoj ) . 
It turns out that for fixed a, the proof is essentially identical to that of Theorem A.l in [DDI) , 
which we briefly reproduce for convenience. For a sequence 7 G F^, we put inil) = li ■ ■ ■ In- 
We recall that if g G PSL(2, C) has large norm, there exists two balls A{g) and R{g) (the 
attracting and repelling balls respectively) on P^, of diameter about jj^, such that g{R{gY) = 
A{g). We also recall that for g G PSL(2,C), 

\ 

-\og\t?{g)-A\ =\og\\g\\ + log 6{g)+0 {I), 

where 5{g) is the distance between the fixed points of g. Thus to control the trace from the 
norm it is enough to control the distance between fixed points. 

Fix e > 0, small with respect to x- Fix a G F, with ||a|| < |n. We have to estimate the 
probability that for large n ^ log tr^(/3(^„(7)a)) ^ [x^ — £,x^l + By the large deviations 
theorem for the distribution of the norms |BLl §V.6] we infer that outside a set of exponentially 
small probability, | ^ log ||/o(^n(7))|| ~Xm| < |- Since ||a|| < ez", we get that with high 
probability. 



n 



log||p(4(7)a)||-XM 



< 



3e 



Therefore we need to show that the probability that (5(p(^„(7)a)) < e 4" decreases like 

g-c(£)n 

For this, let m = [(1 - e)n\ . We split ln{l)a as ln{l)a = (7n • • • 7m+i)(7m • • • lio)- By the 
large deviations theorem for the norms, and the facts that ||a|| < ea"" and m = [(1 — e)nj , we 
get that 

P(|l log \Wm{l)a)\\ - Xm| > e) + P(log Mln ■ ■ ■ 7m+i)ll > 2x^{n - m)) < e-<'>. 

By independence, we may assume that 71 , . . . ,7^ are fixed and satisfy j^log ||/o(^rrt(7)a)|| — Xn\ < 
£ and estimate the probability that (5(/>(£„(7)a)) < e~4" among the 7m+i; • • • ,7n satisfying 
log ||/o(7„ . . . 7m+i)|| < 2x^(ra — m) ~ 2xne. Put h = p{£m{'y)a) and consider the balls A{h) 
and R{h). If we let 



Aip{in{'y)a) = p{jn ■ ■ ■ Jm+i)A{h) and i?(p(4(7)«) 
then by definition we have that 

p{£nh)a) (R{p{inh)ar)=A{p{£nil)a) 



R{h) 



and the diameters of A{p{£n{l)(i) and R{p{in{'y)a) are not greater than e~^'^'^~^'^^'^^'^'^ . So 
if we can show that these balls are separated by a distance > e"*", we then conclude that 
p{£n{'y)a) is loxodromic with a fixed point in each ball, so the distance between its fixed points 
is > e~4" and we are done. But if this does not happen, this means that /o(7„ • • • 7m+i) maps 
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the ball A(h) very close (i.e. at distance < e~4") to the ball R{h). Since h is fixed here, 
the probability of this event is governed by the properties of the unique stationary measure 
on P^, relative to the action of (yo(r),//). It follows that this probability is of order e~^^^ for 
some 6 = 0{p) > (see |DDH Lemma A. 3] for details) and the desired result follows. 

So far we have shown that there exists ci > (which we may assume is smaller than j) 
such that if a is a fixed element in F such that ||a|| < e*", then 



-log|trp(4(7)a)| 
n 



-cien 



Since the number of elements of norm not greater than e^" is equivalent to e'^", we get that 
for C2 < ci, 



3a G r, s.t. ||a|| < e 



C2£n 



and 



-log|trp(^„(7)a)| - Xf^ 



n 



>£ \ <e 



(c2-ci)n 



By choosing C2 = ^, we conclude that p!7j) holds with c 



2 



□ 



Proof of Theorem \3.^ Let us first assume that T is torsion- free. Prom now on is the dis- 
cretization measure for the Brownian motion on X, as defined in ^2.21 We let S = Sr be the 
circle centered at and of radius R and work in the probability space of paths issued from S 
in H, endowed with the probability measure Ws = Wxds{x). For a generic such Brownian 
path, we let T = T{uj) be the stopping time defining fi, Tn = Tn{oj) be the sequence of times 
obtained by repeating the stopping process, and r be the average stopping time, which is the 
almost sure limit r = lim ^ (see the proof of Theorem 12. 7p . Observe that it is enough to 
prove the theorem for the sequence of radii r„ = 

The first step is the following large deviations estimate, which will be proven afterwards. 

Lemma 3.4. Let Tn{u}) be as above. Then for large enough n we have that 



Ws 



n 



> e 



< e 



-c{e)n 



For fixed n, we now define a good set G„ of paths starting from S by specifying its comple- 
ment Bn- Declare that the elements 7 E F appearing in (|17p are bad. From the discretization 
procedure, this corresponds to a set of paths of exponentially small measure, relative to Ws, 
which we denote by B^. Now by definition oj belongs to B^ if \Tn{io) — nr| > where c 
is as in Proposition 13.31 and K is as in Lemma 12.111 By Lemma 12.111 and the strong Markov 
property if w ^ B^, the probability that (iH(w(T„), cj(nr)) is larger than ^ne is exponentially 
small. We let B^ be the corresponding set of paths. Finally we put Bn = B^ U B^ U B^, and 
Gn = B^. Observe that Ws{Bn) < e~<^> . 



We claim that if u is good, that \s u) ^ Gr, 
that 



then for every 7 G Fni^eC 



uj[Tn) 



-ne) we have 



-log|tr(p(7))| - Xti{p) 



n 



< TE, i.e. 



1 



2rra 



log|tr2(p(7))| 

XBrown (P) 



< e 



(recall from the proof of Theorem 12 . 71 that Xfi{p) = ''"XBrown(p))- Indeed by Proposition 13. 3|, if 
oj ^ Bn, then the element q G F corresponding to uj(Tn) (that is, such that uj{Tn) G ol{S)) has 
the property that for every j3 G -Br(«) erne), |^ log |tr(p(/3))| — x^J.{p)\ < t^- Now, if o; G Gn, 
ojirn) lies in -Br(a, %-ne), whence the result. 
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Let E be the set of group elements 7 G Br (^) such that 

> e. 



log |tr2(p(7))| 

XBrown (P) 



2rn 

To complete the proof of the theorem, we have to show that < e'^^^^^^^B^ {^) ■ For this, 
let Gn be the image of G„ under Brownian motion at time nr, that is G„ = {a;(nr), lo E 
The next lemma follows from classical estimates for the heat kernel on M (see below for the 
proof). 

Lemma 3.5. The hyperbolic area of {GnYn -Be (O, ^) is exponentially small with respect to 
the area of B^ (O, that is, there exists c(e) > such that for large n, 

Area ((G„)' n Bm (o, ^)) < exp - c(e)n) . 

We have shown above that if 7 G then -Bh(7, ^1^) is disjoint from Gn- Let 5 = 
where c(e) is as in the previous lemma. By the Vitali covering lemma, there exists a covering 
of £^ by a finite collection of balls {BB_{'yi, 5n),i = 1, . . . A^} such that the balls B^{'^i, |n) are 
disjoint. Thus we infer that 

(18) ^ Area {b^ (^u ^n^) = 47rArsinh2 (^^r^ < Area ((G„)^) . 

Now, a well known theorem due to Margulis [M] (see also [EM]) asserts that when r — )• 00, 

(19) #Sr(r) ^^^^(^=(0'^)) 



Area(A) Area(A) ' 

It follows that there is a constant C such that for every r, ^B-p^r) < Ce^ , hence by (jlSp we 
conclude that 

^ ~ sinh2 {^n) 

and the proof of the theorem is complete in the torsion free case. □ 
Proof of Lemma \3.4\ The almost sure convergence ^ — ^ t was established in the proof of 



Theorem 12.71 the point here is to make this convergence more precise. As we did there we 
slightly shift our point of view and see our Brownian paths as sitting on X. We want to 
formalize the idea of a Brownian path issued from S as obtained from the concatenation of 
path segments starting and ending on S. For this we let fig be the set of continuous paths 
[0,00) — )■ X starting on S, endowed as usual with the probability measure Ws- We define 
a Markov chain with values in Qs as follows: the initial distribution is Ws, and for a.e. cj 
the transition kernel is given by p(w, •) = Waj{T{w)); where T(uj) is our stopping time. This 
Markov chain yields a sequence of random paths (a;„) in rig. By the strong Markov property 
of Brownian motion, the law of the random path obtained by concatenating the path segments 
w„|[o,T(cj„)] is Ws- Associated to the Markov chain there is a real random variable T{uJn) and 
we see that the law of Tn{u) is the same as that of T{ujk)- 

It is easy to see that the random variables (ujn) have the same law Ws- However they are not 
independent since Un+i must start at Un(T{ujn))- The main observation is that independence 
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holds after two iterations. Indeed for every x £ S, 

p{W.r) = I p{a,-)W,{da) = J W^^T(a))W,{da) = J Wyds{y) = Ws 

because by construction the law of a{T(a)) is the circle measure ds. Thus we infer that for 
a.e. uj G S, 



a) = Ws- 



It follows that the random variables a;i,a;3,a;5, . . . (resp. ^2,^4, • • • ) are independent. Fur- 
thermore, by Proposition I2.15"t there exists s > such that E(e'^^('^'') < 00. Therefore the 
desired estimate follows from the classical large deviation theorem for partial sums of i.i.d. 
random variables (applied separately to X^even k'^i^k) and Ylodd k'^i^k))- D 

Proof of Lemma \3.5\ Let vt (resp. ff) be the law at time t of Brownian motion with initial 
distribution ds (resp. 5q). For notational ease, we put F = n (0, ^) . We know that 
^nriF) < Ws{Bn) < Q-c(e)n ^ from this we want to derive an estimate for the hyperbolic 
area of F. Let J > to be specified later, and decompose F as 

(TIT' \ \ / / TIT TIT \ \ 

0, — -Ti6)) U i^FnAM[—-n6,—)) 

where Aji{r, s) denotes the annular region centered at with inner radius r and outer radius 
s. The first observation is that 



n5] . 



Area (f n Bm (o, ^ - nj) ) < Area (o, ^ - nd) ) < exp 

The probability measure is radial, and from the estimate for its radial density k^{t,r) 
given in ([U , we deduce that for t = tit and ^ — nS < r < ^ , 

U-"i^<A;0(nr,r)<^. 
n X n 



Since v^r is an average of bounded translates of f^^, the same estimate holds for the radial 
density k{nT,r) of i^nr- On the other hand the radial density of hyperbolic area is rsinhr, so 
we infer that 

(/ TIT TIT W {' 

F n Am ( tl5^ — ) ) < / r(sinh r)(ir(i0 

^2 2 / / J FnA^{^-nS,^) 

< / r (sinh. r)k{nT, r) \/ne"'~ drdO 

< T?!'^ exp + ^"-^ VnriF) < n^/^ exp + - c{e)n 

Finally we get that 

Area(F) < exp - nb^ + n^^^ exp + - c{e)Ti \ , 



so by choosing b = ^/rde) we are done. □ 
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Let us now assume that F has torsion, and briefly explain how to adapt the proof of the 
theorem. We pick a torsion-free subgroup F' < F of finite index, and apply the first part of 
the proof to F'. In this way we get a discretization measure /i on F'. We slightly modify 
Proposition 13.31 to allow for a perturbation a belonging to F (this only slightly affects the 
counting in the last part of the proof), so that its conclusion now reads 



(J^-f G F', 3a G F s.t. ||a|| < 



and 



1 



n 



log|trp(7a)| - Xfi 



> e > I < exp(— cen). 



Then we argue exactly as in the torsion free case, by introducing the set of 7 G Fni?H (^) 
such that I 27;^ log |tr^(yo(7))| — XBrown(p)| > To the discretized Brownian motion on F', we 
associate the sets Gn and G„ exactly as before, and we estimate #E by arguing that, from 
the modified version of Proposition 13.31 if 7 G then -Bh(7; ^n-) is disjoint from G„. This 
completes the proof of the theorem. □ 
The following corollary, of independent interest, is the key step in the proof of Theorem 
311 



Corollary 3.6. Let T be a lattice in PSL(2,M), and p be a non- elementary type preserving 
representation ofT. If {"fn) is an admissible random sequence in T, then almost surely, 

— ^log|tr2(/9(7„))| > XBrown(/3)- 

Proof. Recall that admissible means that the 7„ are chosen randomly in Br{rn), relative to 
the counting measure, where r„ is a sequence such that ^ e"'^'"" converges for every c > 0. 
Fix e > 0. We will show that 



(20) 



1 



#i?r(r. 



# 7 G Brir„ 



1 



4d(7„) 



log|tr2(/)(7„))| 

XBrown 



> e 



We then deduce from the Borel-Cantelli lemma that a.s. for large n, ^^^^ ^ log |tr^(/3(7„))| is 
e-close to XBrown (p) < and the result follows. 

The first observation is that in BY{rn), by the Margulis estimate (|19p on ^B^ir), for every 
6 > 0, we have that 



(21) 

We choose 6 = min ( ^ > ^ 



#{7, d(7n)<(l-'5)r„}<e(i-^>". 
. Now if 



diin) > (1 - S)rn and 



■log |tr2(p(7„))| 

XBrown (P) 



an elementary computation shows that 



1 



^log|tr^(p(7„))| - XBrown (p) 



> (1 - 6)e - (^XBrown (/>) > ^• 



Therefore the desired estimate (|20p follows from (|21|) . Theorem 13 . 2 1 and the fact that ^ B^irn) 
grows like e^". □ 

Remark 3.7. This corollary yields yet another definition of the natural Lyapunov exponent, 
where the choice of Riemann surface structure is reflected in the fact of counting in the 
lattice F. Similarly, we could replace the trace by the norm, and assert that generically. 
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2^^^ log 11/3 (7„) II converges to XBrown(/5)- The proof is identical to that of Theorem 13.21 
(observe that Proposition 13.31 becomes much simpler in this case) . 

3.2. Proof of Theorem 13.11 Let rrin be the normalized counting measure on BY{rn) and 
m = nn>i™'"- considering the potentials of the currents ^^^^[^(7^, i)], it is enough 

to show that for m-a.e. sequence 7 = (7„) in T^, the sequence of psh functions (A 1— 
M(A,7,n))„>i on A defined by 

u(A,7,n) = ] log \i?{pxirin)) - t\ 

converges in L\^^{K) to A 1— )• XBvown{p\)- Corollary 13.61 above implies that for every fixed 
parameter A, u(A,7,n) converges m-almost surely to XBrown(PA)- Now clearly the set E C 
A X of parameters (A, 7) such that ti(A,7,n) converges to XBrovmiPx) is measurable, so 
by Fubini's theorem it has full d\®m measure. Reverting the order of integration (Fubini's 
theorem again), we infer tlia-t for vn a,.e. 'ii(A, 'y, tz) converges to XBrown {px) dX-a.s. 

In addition, the sequence n(A,7,n) is locally uniformly bounded from above on A. Indeed 
we have that 

log|tr2(p;,(7„))-t| < ^^log||/,;,(7„))f + 0(1) 



Mln) " ' " ' ' - M7„) 

^2!^)^°^"^""+^^'^-^+''^'^ 

where the constant f3p^ is locally uniform in A (see Corollarv l2.3p . By the standard compact- 
ness theorems for sequences of subharmonic functions (see Theorems 3.2.12, 3.2.13 and 3.4.14 
in [Ho] ), we conclude that n(A,7,n) converges to XBrown(pA) in -^ioc(^)- Indeed, if 7 is such 
that convergence holds a.e in A, we can extract a subsequence converging to some v in Lj^^, 
and limsup„_^(^ u(A, 7, n) = v outside a polar set, therefore v = XBrown- This finishes the 
proof. □ 

Remark 3.8. As for Theorem 13.21 the statement of Theorem 13. II mav be adapted to the case 
where T admits non-trivial torsion. 

3.3. Random closed geodesies. Let CG{X) be the set of closed geodesies in X, which is 
a countable set. We allow a geodesic to be travelled several times. If this does not happen, 
the geodesic is said to be primitive. 

Let p G Hompar(G, PSL(2, C)). If 7 G CG(X) the holonomy hp('y) is defined up to con- 
jugacy in p{T), so it makes sense to talk about its (squared) trace tr^/ip(7). If now [px) is 
family of parabolic representations, for t G C we may consider the subvariety in A defined by 

Z(7,t) = {AG A, tvHp,{j) = t} . 

Here we first briefly review several notions of "random closed geodesies" on X, and interpret 
Theorem 13.11 in these terms. 

3.3.1. Thurston model. To start with, assume that X is compact. The following model for 
random geodesies on X was apparently popularized by Thurston (see e.g. [BoJ): for {x,v) G 
SX let as before "y(x,v) be the associated unit speed geodesic ray. For t > 0, close the path 
7(x,t;) I [o,i] with a shortest path returning to x. Notice that there may be a choice involved here, 
which we may disregard for it happens only for a set of {x, v) of zero measure. We denote by 
cg(x, v) (t) the unique closed geodesic freely homotopic to the obtained loop. Notice that its 
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length may be much smaller than t. Nevertheless it follows from Theorem 13.21 applied to the 
identity representation that this phenomenon happens only on a set of exponentially small 
measure. 

In this way we have constructed a family of projections v i— )• cg(x, v){t) from SX to CG{X) 
indexed by t. Projecting the (normalized) Liouville measure gives rise to a family of proba- 
bility measures (?Tif°°*^''^"^)f>o on CG{X). 

Several variants of this construction may be considered: 

- Variant 1: a first possibility is to restrict to geodesic paths starting from the marked 
point -k, and close them so as to obtain loops based at 

- Variant 2: allow for paths starting from any {x,v), but append them at the two 
endpoints with shortest paths joining them to 

Notice that in these two variants, we naturally obtain a family of measures on the fundamental 
group, which in turn projects to a family of measures on CG{X). Notice also that in the 
second variant, the family of measures that we obtain in CG{X) is the same as in the original 
Thurston construction (since the corresponding loops are conjugate in T). In any case, the 
measures {'m^'^"'^'^^^'^)t>o may be thought of as measures on F. A good intuitive picture for 
^geodesic ^^^^ ^£ being loosely equidistributed on an annulus of bounded width about the 
circle of radius t in F C H. 

These variants are interesting because they also make sense when X is only of finite volume. 
Indeed, using these constructions we obtain a family of measures (?7if*^°'^'^^'"^)t>o on G. Now 
7 G F projects onto a closed geodesic if and only if it is not parabolic. Again by Theorem 
13.21 this happens outside a set of exponentially small measure. So what we do is simply to 
restrict to loxodromic elements, renormalize the mass and project to CG(X) we also obtain 
in this way a reasonable model for random closed geodesies on X. 

The following lemma asserts that the measures (mf*'°'^'^'^"^) are well spread in F. 

Lemma 3.9. Let X be a hyperbolic Riemann surface of finite type endowed with its hyperbolic 
metric, and let {mf'^°'^'^^^^)f:>o be any of the family of measures on F constructed above. Then 
for every 6 > there exist constants K and c such that for all t > 0, for every subset A (ZT 
with cardinality < e^^"'')*, we have that mf'^°'^'^^^'^ (A) < Ke~'^* . 

Proof. Let us treat the case of paths starting from -k (Variant 1 above). If is a unit tangent 
vector at we denote by 7„ the unit speed geodesic ray starting from * in the direction of v. 

The first claim is that there exists a constant Ki depending only on X such that the 
probability that dxilvit),-*^) is larger than |t is not greater than K\e~^^l'^ . A proof of this 
fact is outlined in [Sul', §9], and discussed with much greater detail and generality in [KMj . 

For t > 0, let Et be the set of tangent vectors v G T\X such that dx{riv^\'^) ^ f*! which 
by our previous claim has measure > 1 — K\e^*'l'^ . We lift the situation to H, so that our 
geodesies are now rays of hyperbolic length t issued from 0. If a lattice point 7 G F • is the 
endpoint of 7„(t) for some v G Sj, then c?H(75 7t;(0) < We apply the first cosine rule [Beat 
§7.12] in the triangle T(0, 7, 7i;(t)) to get an estimate of the angle d at the origin between the 
rays 7^(t) and [0,7]. It follows that cos(6l) = 1 - Ce-(2+°)* + 0(6"^*), with |a| < (5, from 
which we deduce that d < K2e~^^~^^'^^^ . Hence for any subset A of F, the set oi v £ Et such 
that 7^(t) belongs to A has measure at most i^2(#^)e~^^~^/^^*. 
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We conclude that if A is as in the statement of the lemma, the measure of the set of 
tangent vectors v such that ^v{t) belongs to A is not greater than [Ki + i^2)e~2* and the 
result follows. □ 

It is an entertaining exercise in hyperbolic geometry to show that the cardinality of the 
set of 7 E -Br(t) which are proper powers of elements of T is of order of magnitude e*/^. It 
follows that a generic geodesic with respect to the measure TTiSeodesic primitive. 

3.3.2. Brownian Thurston model. We can do exactly the same construction by using Brownian 
paths instead of geodesic ones. As before, several variants are possible, by choosing the initial 
point to be either fixed or distributed according to the Poincare volume element on X, and 
by choosing the paths to be closed to become a loop based at * or not. In this way we 
obtain families of measures on F or on CG{X), which, abusing slightly, will be given the 
same notation (772-|^'^°™'^'^'^)t>o- In the unit disk picture, the corresponding measures on T are 
induced by the heat flow on H^, so they now look like Gaussian measures concentrated around 
the circle of radius ^t. 

Again, these measures are well spread on T. 

Lemma 3.10. Let X be as above and {mf^°'"^^^^)tyo be the family of measures on T obtained 
from the Brownian model for random geodesies. Then for every 6 > there exist constants K 
and c such that for all t > 0, for every subset A C T with cardinality ^A < e^^^^^^^ it holds 
that mp°^'^'^'^(yl) < i^e~^*. 

The proof uses arguments similar to those of Theorem 13. 2|, and will only be sketched. 
Details will be left as an exercise to the reader. 
We need the following estimate. 

Lemma 3.11. Let {X,-k) be as above. There exists a positive constant c such that 

P,(dx(*,^(t)) >ct)<e-^*/2^ 

Proof (sketch). Let / be the density of the diffusion of 5^ at time 1, relative to the hyperbolic 
area element on X. The classical Li-Yau estimates on the heat kernel (see e.g. |CC1 §B.7]) 
imply that / € L°°{X) C L'^{X). Let {Ilt)t>o be the Brownian semigroup on X, and put 
Vt = {x & X, dx{x,-k) > ct}. We have to show that Jy^Htf is exponentially small with t. 
For this, using the Schwarz inequality we write 



Vt 



<l|nj||i2(x)Vol(T4)'/' 



and we conclude using the fact that ||IIi/||^2(x) converges, while vol(Vj) decreases exponen- 
tially with t. □ 

Proof of Lemma \3.10\ (sketch). We lift the situation to H, and look at Brownian paths issued 
from 0. Fix a constant c > 0. By the estimates ^ for the heat kernel on EI and the above 
lemma, the probability that io{t) ^ A^ {^{1 — c), ^{1 + c)) or dM{(^{t),T) > ct decreases 
exponentially with t. We discard this set of paths, so that it is enough to work with A' = 
ACiA^ (1(1 — 3c), 1(1 -|- 3c)) instead of A. Now if we denote by uit the concatenation of u;|[o^t] 
with a shortest path joining uj{t) to F, we infer that P {ujt £ A') < vt (UaeA' -^(^j ci)) , where 
Vt is the law of Brownian motion at time t. We compare this probability to the hyperbolic 
area of IJaeA ^{o., ct) by using ^ again, and get an estimate of the form P (cot G A') < 
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exp((— (5 + 5c)t). Finally, choosing the constant c small enough with respect to S finishes the 
proof. □ 

3.3.3. Length-based model. Let CGt(X) be the set of primitive geodesies of length at most t 
on X, and let 171^"^^^^ be simply the normalized counting measure on CGt{X). 

To apply Theorem 13.11 we need to to lift rn^'^^^ to F. For this, we choose a shortest path 
from ★ to a given geodesic to make it a loop based at Observe that the length of this path 
is uniformly bounded. Indeed, recall from Proposition 12.21 that there are neighborhoods Ci of 
the cusps such that each geodesic must enter the compact region X\[jCi. 

We then infer that the lifted measure is supported in an annulus of bounded width about 
the circle of radius t in F. Also, the lifted measure is well spread in F in the sense of Lemmas 
13.91 and 13.101 because of the well-known asymptotic estimate 7^CGt(X) ~ y as t — )• oo. This 
estimate originates in the works of Huber and Selberg (references covering the non-compact 
case include |Hel H]) 



3.3.4. Equidistribution. We now choose any of the above models for random geodesies, and 
denote by (m^) the associated family of measures on CG{X). Let us recall the statement of 
the equidistribution theorem on A associated to random closed geodesies (Theorem [U]) . 

Theorem 3.12. Let X be a hyperbolic Riemann surface of finite type, and {p\)xeA 
holomorphic family of representations of G = iTi{X), satisfying (Rl, R2, R3). Let Tbif = 
XBrown ^6 ^hc natural bifurcation current. 

Let (tn) be a sequence of positive numbers such that for every c > the series Yl e^'^*" 
converges. Let (7^) be a random sequence of closed geodesies on X , according to the product 
measure W rrit^ . Then almost surely we have that 

l—-[Z{^n,t)] Tbif. 

41ength(7„) " n~^oo 

The proof of the theorem is identical to that of Theorem 13.11 using the fact that the lifted 
measures mt„ are well spread on F. Indeed for fixed A, a Borel-Cantelli argument based on 
Theorem 13.21 together with Lemma 13.91 (Thurston model). Lemma 13.101 (Brownian Thurston 
model) or the estimate ^CGt{X) ~ ^ (length-based model) implies that almost surely 

— 1— ^log|tr2(/9(7„))| > XBrown (/?), 

41ength(7„) ' ' n-^oo 

which is the statement of Theorem |Dj Theorem 13.121 (i.e. Theorem [C]) then follows exactly 
as in ^3.21 



4. Further remarks, open problems 

4.1. A first natural problem is to study the dependence of the natural Lyapunov exponent 
XBrown ^s a function of the Riemann surface structure X. It is likely that this function is 
continuous in {p,X). It would be interesting to study the possible subharmonicity properties 
of this function as X varies in Teichmiiller space. 
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4.2. (cf. jPDH §5.2.3]). The varieties Z('y,4:) appearing in the equidistribution theorems 
are made of accidental parabohcs and new relations. The question is to determine if one 
of the two types dominates. At first sight one might guess that the varieties of accidental 
relations {X, pxi'j) = id} are of codimension 3 in the representation variety. Unfortunately 
this is wrong. Indeed, if tr^(p(7)) = (a codimension 1 condition) then p(7^) = id. More 
generally for any k we can obtain hypersurfaces of representations such that p{'y^) = id. 
On the other hand, as observed in H'6.'6\ random geodesies are generically primitive, so this 
phenomenon is negligible (see |Luj for a wide generalization of this fact). 

Is it true that for a generic sequence (7„) of closed geodesies, {X, px{['yn]) = id} is asymp- 
totically of codimension greater than 1 in the representation variety? 

4.3. Our results may be recast in the more general framework of holomorphic families of 
PSL(2, C)-cocycles over a discrete or continuous time dynamical system, endowed with some 
invariant measure. A general question in this area of research is how the values of this cocycle 
over the set of periodic orbits approximates the Lyapunov exponent (see e.g. [Kal] ) . What 
Theorem [D] says is that we can have very precise such results (with exponential deviation 
estimates) in the particular case of the geodesic flow in constant negative curvature. It 
would be interesting to know how this can be generalized to a more general (say, uniformly 
hyperbolic) base dynamics. One possible way of doing this would be to use the theory of 
Markov driven random products of matrices, which was developed in particular by Guivarc'h 
and his collaborators (see |Gu] for a review). 

4.4. What can be said about the distribution of [Z(7„,4)] when the 7^ are simple closed 
geodesies on X7 
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